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Abstract 



We use dynamical systems theory to construct the normal form 
^ • of the Navier-Stokes equations for the flow of a thin layer of fluid 
t> upon a solid substrate. The normal form equations illuminate the 
^ fluid dynamics by decoupling the long-term flow of interest from the 
^ ' rapid viscous decay of the transient shear modes. The normal form 
clearly shows the centre manifold of the lubrication model and shows 
the result that the initial condition for the fluid thickness of the lu- 
brication model is not the initial physical fluid thickness, but instead 
is modified by the initial lateral shear flow. With these initial condi- 
tions, better forecasts will be made using the lubrication model. This 
dynamical systems approach will also enable similar illumination of 
other complicated models of dynamics. 
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1 Introduction 

The fiow of a thin layer of fiuid is important in numerous industrial and 
natural processes: industrial applications include the coating processes of 
autobodies, beverage containers, sheet goods and films, decorative coating, 
and gravure roll coating; biomedical applications include the liquid films 
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covering the cornea of the eye or protecting the hnings of the lungs A 
variety of mathematical models of the fluid dynamics have been developed to 
aid understanding and prediction of the flow. For examples see the extensive 
reviews of Chang [2J, Ruschak 15^ or Oron, Davis and Bankoff On a flat 
substrate and in one lateral dimension the usual model for the evolution of 
a viscous fluid film's thickness ri{x,t), driven only by surface tension, is the 
following nondimensional "lubrication" model |2l e.g.] 



where the Weber number We measures the strength of surface tension in the 
thin film flow. This is a considerable simplification of the governing equations 
when compared to the full Navier-Stokes equations. 

Consider making a forecast with the lubrication model from some 
initial state of the fluid, say with physical thickness rjo. Surely one just sets 
the initial field ri{x, 0) to the initial thickness of the actual fluid t]q{x). Herein 
we discover this is incorrect and that instead the correct initial condition for 
the lubrication model ((T)) is more accurately 



where uo{x,y) is the initial lateral velocity field of the fluid, and Re is an 
appropriate Reynolds number. This initial condition accounts for the flow of 
fluid in the transient viscous decay of the initial velocity field. 

For example, two colliding sheets of fluid form a hump. Suppose the 
initial fluid state is that of a constant thickness, ?7o constant, and Couette 
flow, Mo = Uof{x)y/r]Q where f{x) describes the lateral variations in strength 
of the initial lateral shear flow. Then from (jSJ the corresponding initial 
thickness for the model is 




(1) 




(2) 



?7(x, 0) ~ ?7o — Re 



24 



(3) 
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Two colliding sheets of fluids corresponds to a localised and negative deriva- 
tive f'{x); in this case the initial condition Q says the best forecast is 
obtained by starting the model from an initial condition with a hump pro- 
portional to — /'(x). This hump accounts for the transient collection of fluid 
at the collision as the lateral velocity decays through viscosity. This physi- 
cally reasonable result agrees with earlier work ^3 • 

Here we develop analysis, Section 13. 3| to determine that the parabolic 
weight function is correct in the integral of Further, we introduce and 
use a framework. Section 13 to show that Q is just the flrst approximation 
in an asymptotic series for the correct prescription of the initial conditions 
for the lubrication model (Q), see Section 0J The theoretical framework is 
that of the normal form transformation in dynamical systems with centre 
manifolds jH El . The continuity equation of incompressible fluid dynamics 
is an algebraic equation; thus, we make the novel application of the normal 
form transformation to differential algebraic systems. 

Because of the complicated detail and form of the fluid dynamics equa- 
tions, we introduce the normal form transformation in Section El using a 
relatively simple example dynamical system. The transformation is a near 
identity transformation with the aim of nonlinearly decoupling the dynami- 
cal equations as far as possible. In systems with a centre manifold, the centre 
manifold model appears immediately, see the example of Section 12.11 as flrst 
recognised by Elphick et al. Cox and Roberts then recognised that 
the normal form also immediately provides the correct initial conditions for 
a model, see the example of Section IT^ The example is a coupled differen- 
tial algebraic system; in Sections 12. 3H2 .41 we show that the normal form may 
none-the-less be constructed in the usual manner. 

Then in Section |21 we introduce the fluid equations and develop the nor- 
mal form transformation that decouples the viscously decaying dynamics of 
lateral shear from the long lasting slow evolution of the thickness of the fluid 
fllm. When decoupled, the arguments of Section |3] provide initial conditions 
for the lubrication model (^. Equation Q gives the leading nontrivial terms 
in the asymptotic expression for the initial condition. 
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2 An example introduces the normal form 

Here we explore a straightforward artificial example (introduced in [16J) to 
illustrate the techniques used later in the normal form derivation of the initial 
conditions Q for the lubrication model (Q). This example reflects in detail 
many essential aspects of the relation between the normal form and initial 
conditions, and how to use the freedom inherent in the derivation. Someone 
very familiar with normal form transformations may skip this section. 

Consider the following index 2 differential-algebraic system which is struc- 
turally analogous to the thin fluid film equations (jlHHSDI) introduced in Sec- 
tion ini (variable subscripts denote differentiation) : 



Some trajectories of solutions are shown in Figure ^ In this section u{t), 
v(t), pit) and rjit) represent analogues of a fluid film's velocity, pressure 
and thickness. The small parameter e is introduced to mimic slow lateral 
variations in the fluid thickness and flow. The normal form transformation 
we explore in this section clarifies the long term evolution of the above finite 
set of ordinary differential evolution equations. 

Low-dimensional dynamics: When e = , the system (jlHZI) has the equi- 
librium solutions rj = const, u = v = p = (mimicking a uniform motionless 
fluid layer of thickness 77). These equilibria, call the set of them A4o, are sta- 
ble: linearising, putting e = and seeking nontrivial solutions proportional 
to exp(At) leads to requiring the determinant 



Ut + ew^ + uv + ep + (l + e^) m + e?7 = , 
Vt + euv + + p + (1 + e^) f = , 



(4) 
(5) 
(6) 
(7) 



ew + f = , 
Tjt + euTj — f = . 



A + 1 








A + 1 
1 

-1 





1 







A 



A(A + 1) = 0. 



(8) 
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Figure 1: simulations of the example system (@HZI) for e = 0.4 over time 
interval 6 at time steps of At = 0.2 and projected onto the r/u-plane. This 
shows the exponential quick collapse of states onto the low- dimensional centre 
manifold (solid curve) given by ©. 



This shows: the rapid (9(e~*) decay to any of the equilibria is determined 
by the u-equation (jlj); equation ((Tj) determines r] is a. neutral mode; whereas 
the other two equations do not contribute dynamics because of the need to 
satisfy the algebraic 'continuity' equation (jHl). For small e 7^ , instead of 
7] being constant, we thus expect that 77 will evolve slowly on a nearby and 
similarly exponentially attractive set of states, called the centre manifold. 
This is confirmed by the numerical simulations shown in Figure ^ We soon. 
Section 12.11 confirm this nearby centre manifold of slow long term evolution 
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IS 



u 



p 



V 



er^(l-3e2)+eV + (^(e'), 



(9) 
(10) 

(11) 



On this centre manifold the slow model evolution is 



r]t = e\l + r])r] - e\l +r]){3 + r]) + O(e^) . 



(12) 



By the centre manifold relevance theorem 7,, Chapt. 5, e.g.] or [Ij, we expect 
the long term behaviour of all nearby solutions are described by (jHFEI)- 

Elphick et al. |4| noted that low-dimensional centre manifold models such 
as (jnHT^ are an immediate consequence of transforming differential equations 
into a normal form. Cox & Roberts [3j observed that the projection of initial 
conditions onto such models also immediately follows from such a normal 
form. However, neither group explicitly addressed dynamics governed by 
systems of differential-algebraic equations. Here we demonstrate how the 
techniques adapt easily to differential-algebraic systems such as both (j^HZI 
and also incompressible fluid dynamics, see Section |S1 We illuminate the 
dynamics and modelling of such systems with a normal form transformation. 

2.1 Interpret the normal form transformation 

In Section 1231 we will show how to construct the near identity transform to 
new dynamical variables that simplifies the description of the dynamics. But 
first we show how the resultant transformation illuminates the dynamics by 
decoupling the long-term evolution from the short-term decaying transients, 
and by decoupling the algebraic component of the governing equations. Sub- 
stitute into the dynamical system (pMl2j). with new variables denoted by 
Fraktur font, the near identity transformation: 



V 



+ e(l + t))u + e'Kl + f))u2 + e=^(l + f)) [-(2 + [))u + ^ ] 
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+ e\l + m-il + IhW + >1 + 0{e') , (13) 
u = u - ef) + e'[{3 + f))f) + i(l + [j^] + e%{l + f))u=^ + 0{e') , (14) 

V = 0-eu + e2(^-e^[(3 + [))[) + i(l + [))u2] + O(e5), (15) 

p = p-e%-e%l + l))u + e%-l{l + [)y] + 0{e'^). (16) 
Then the example dynamical system (jlHIj) becomes 

[), = 6^(1 + f))f) - e^(l + f))(3 + f)) + 0(e5) , (17) 

= - [ 1 + e\2 + f)) - e^(l + f))(3 + f^) ] u + O(e^) , (18) 

D = p = 0. (19) 

Crucially, the normal form system l{r/\\iy\} captures all the solutions of the 
original system at least in some neighbourhood of the equilibria A4q. 

The reason is simply that the transformation p3flTH|l is a smooth reparametri- 
sation of the complete state space near Aio. Thus from arbitrary feasible 
initial conditions the normal form system (jl7H19|) retains all the transient 
dynamics and all the long-term dynamics provided the dynamics stay within 
the neighbourhood of A4o- Consider the effects of this transformation. 

1. The algebraic part of the original system (0HZI) is transformed into the 
trivial and decoupled algebraic equations (fTIH) , namely = p = . 

2. See immediately from the form (|18p that u = is an invariant manifold 
of the dynamics, and is exponentially quickly attractive to a wide vari- 
ety of nearby initial conditions. Thus this normal form clearly displays 
the attractive centre manifold is u = . In original variables, the trans- 
formation (jl3H16|) then immediately shows the centre manifold maps 
to ?7 = f) and the earlier claimed (jUHTTj). 

3. Also see the evolution (jl7j) of i) reaffirms the earlier claimed centre 
manifold model (jl!2j) . 

4. Lastly, a crucial feature of the {) evolution (fTTj) is that it is independent 
of u in this normal form. Thus all states with the same f) but different u 
only differ in the evolution of u, the f) evolution is identical. Next we 
use this to deal with initial conditions. 



Tony Roberts, February 8, 2008 



2 At 



0.742- 
0.594- 
0.445- 
0.297- 
0.148- 
0.000 
-0.148^ 
-0.297- 
-0.445- 
-0.594^ 



-0.742 




-1.05 -0.84 -0.63 -0.42 -0.21 0.00 0.21 0.42 0.63 0.84 1.05 



Figure 2: centre manifold u = (thick curve) of the example system (@HZI) 
for e = 0.4 . The 16 transverse curves, called isochrons, show initial points 
which have the same long term dynamics. 



2.2 Project initial conditions 

Consider the result of specifying some initial condition in a simulation, say 
v{0) = r]o, u{0)=uo, v{0) = vo, p(0) = po • (20) 

In general these will not lie on the centre manifold model (jnHTT|l. see for 
example the simulations shown in Figure ^ Following the arguments in jS], 
we use the normal form to deduce how to project such an initial condition 
into one for the low- dimensional model ()12|l . 

But first we numerically illustrate the projection. Numerical solutions are 
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feasible for simple systems. By running simulations, such as those in Figured 
we draw contours of all initial conditions whose subsequent solutions end up 
at the same location (to an exponentially small difference) on the centre 
manifold after some fixed time, see Figure |21^ AH the points reaching the 
same location at the same fixed time must then have the same long term 
evolution. Hence, all initial conditions lying on one of the contours in FigureEl 
must be projected onto the centre manifold where the contour intersects the 
manifold. Then the model will predict the correct long term evolution from 
the originally specified initial condition. 

Any such initial condition corresponds to some point (f)o, Uq, Oq, po) in the 
transformed variables. First, revert the series in the transformation (|13H16p 
to find the initial state in the normal form variables 

ho = ^0 - e(l + ^o)mo - e^(l + ^o)(^o - |mo) 
+ 6^(1 + rio) {18 + l8rio-ul)uo 

+ e^^(l + Vo)i8Ar]^o + 144r7o - gOr^o^^ - 102ul + 4) + 0{e') (21) 
Uo = Mo + e?7o-e^(l + ?7o)Mo-e^(4 + 2r7o)r7o 

+ ie^(l + r/o)(24 + 4 + lQr]o)uo + 0{e') , (22) 
Oo = ^^0 + ewo , (23) 
Po = Po + e^Vo-e% + 0{e'). (24) 

Second, for this system with an algebraic component the initial conditions (j2(Jj) 
must be consistent with the algebraic constraints (fTTHl . namely Do = po = ; 
hence, we require the original system to satisfy 

vo = -euo and po = -e^m + ^^Vo + C'(e^) . (25) 

Finally, when Uq is non-zero the initial condition is off the centre manifold 
u = . But in the normal form ()17l[TK|) the evolution of f) is unaffected by 
the evolution of u. Further, u decays to zero exponentially quickly. Thus, 

^ These contours are called isochrons [01121; they are also known as the stable fibra- 
tions "S', §5.2]. 
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apart from an exponentially decaying transient, the evolution from {i)o,Uo) 
will be identical to that from the point (P)o, 0) on the centre manifold. Thus 
the appropriate initial condition for the low-dimensional model (jl2|) is not 
ri{0) = TjQ but instead that 7/(0) = f)o as given by (j2H). This projection of the 
initial condition takes into account the initial transients in the dynamics. 

2.3 Construct the normal form 

We now explore how to construct the particular normal form transform (jl3l - 
ITB|) that has the requisite properties we used in the previous Section IT^ The 
construction has considerable detail, see [SIIIHIE], with many subtleties in 
its application to this problem. To find the transformation (fT!^lfi|l with cor- 
responding evolution ()17flT!?|) . we seek a near identity coordinate transform 
in the general form 

1] = i) + H{[),u,e) , M = u + f/(f),u,e) , 

v = v> + VH),u,e), p = p + P(f3,u,e), (26) 

such that the differential-algebraic system (@HZI) takes this differential-algebraic 
normal form 

f)t = 7^(f),e), Ut = -u + W([),u,e), D = p = , (27) 

where the functions H, U, V, P, Ti and U are strictly nonlinear functions of 
f), u and e, and where we require 

U strictly linear in u , that is W oc u . (28) 

That these properties (j2ZH2HI) can be found is assured by the linearisation 
leading to (jHl): the zero eigenvalue assures that // ~ f) evolves nonlinearly; 
the —1 eigenvalue assures us that m ~ u decays linearly with some nonlinear 
modifications; and the lack of other eigenvalues assures us that the normal 
form has two algebraic constraints on the other variables. The condition (j28|) 
ensures that u = , with = p = , immediately describes the exponentially 
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attractive centre manifold Consequently, the model evolution is simply 
[)t = 7i(f), e) , from (j77|) . and the projection of initial conditions is done along 
curves of constant f). 

The most straightforward construction uses iteration similar to that in- 
troduced for constructing centre manifolds jTTj. Suppose we have some cur- 
rent approximation to the transformation and the evolution; for examples, 
H = V = U = P = T-l = V( = Oisan initial approximation, whereas (|13H18p 
is an approximation with errors (9(e^). Given some such approximation, we 
seek small corrections, denoted by primes, that improve the current approx- 
imation: that is, for example, we seek u = u + U + U' for some small U' to 
be determined. Substitute these into the governing system (@HZI) and omit 
products of small corrections to obtain a system of equations to solve for 
corrections. For example, the left-hand side of the ?7-equation ((Zj) becomes, 
upon using the chain rule, 

rjt + euT] — V 
drj di) dr] du 



1 + 



OH dH'\^^, fdH dH'\ ^ 



+ e{tx + U + U')(lc, + H + H') - + V + V) 

1 + ^ 1 H + ^(-u + U)+ e(u + U)i[, + H)-ix, + V) 

ou 
f)ff' 

where denotes the residual of (jZj) at the current approximation. Hence 
to choose corrections which reduce the residual we need to solve (|32j) below. 
Similarly for the other equations to require the solution of 



u- 



dir_ 

du 



U' = i^+W, (29) 
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(30) 
(31) 
(32) 

where Ri denotes the residual of the system's equations (jlHII) evaluated at 
the current approximation. For example, at the initial approximation H = 
V = U = P = Ti. = U = 0, and using D = p = , the residuals 

= eu^ + e^u + eP) , i^=0, /^=eu, and /^=euf). (33) 

Now, see that the algebraic equation (jHl) immediately gives the ^-correction 
through (p?T|l . and that the dynamic equation © for v immediately gives the 
P-correction through ()3()|1 . For example, in the initial iterate V = — eu and 
P' = . The other two components are more subtle. 

• For equation (|32j). since V is known, any term involving a factor u'' in 
the residual V will cause, via the homo logical operator u^-, a 
term with u'^/q in the transformation H'; for example, with V = 
eu(P) + 1) the correction H' = eu(P) + 1) . However, any term in the 
residual -f^— V with no u factor cannot be matched by a term in H', 
as u^- necessarily involves u, and must instead be matched by a term 
in the evolution W. Thus the transformation is found so that the 
evolution of P) only involves f) itself as required by (^Tjl . 

• For equation (j29j) any term involving a factor in the residual ifeiwill 
cause, via the homological operator — U', a term with uy\q — 
1) in the transformation U'; for example, with = eu^ + el) the 
correction U' = eu^ — ef) . However, any term linear in u, such as the 
term e^u in (fHHjl . must be used to correct the evolution via W. Thus 
the transformation is done so that the evolution of u satisfies (|28|) . 

In this manner iteration builds asymptotic approximations to the normal 
form transformation and the evolution, such as the earlier (fTHMl9|) . 
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2.4 Homogeneous solutions provide freedom 

There are extra complications because the normal form transformation is 
not unique. In general the transformations ()2fij) are chosen to minimise the 
number of terms in the corresponding normal form dynamical equations ()27|l 
by placing as many terms as possible in the transformation ()26p . However, 
the result then depends upon our choice of basis for the algebra. Here we 
briefly explore the options available to us while retaining all the properties 
required to extract the centre manifold model and the projection of initial 
conditions. 

As an aside note that the essential properties of the normal form are 
displayed in Figure |21 the curve of u = must correspond to the centre 
manifold (the thick curve); and the curves of constant f) must correspond to 
the isochrons (thin contours). Freedom comes from being able to parametrise 
these curves in state space in any smooth way consistent with these require- 
ments. 

The simplest way to identify freedom is to find the homogeneous solutions 
of the correction equations ((2nHS21) for H, V, U, P, 7i and U. The full 
homological operator jHl e.g.] for this problem appears on the left-hand side 
of (I2nHS2I); it operates on the space of multinomials in u, f) and e, but only 
the polynomials in u are significant as f) and e do not appear in the operator. 
First see that homogeneous solutions must have V = P' = Q . Then, second, 
two families of linearly independent solutions are 

{H\n') = {v.\qu''), U'=U' = 0, (34) 
and (f/', W) = (u^ (g - l)u«) , H' = H' = . (35) 

Thus we can always change an iterate by some linear combination of these 
solutions, multiplied by some function of f) and e, and only affect higher-order 
terms in the transformation. 

• Consider (j34|) with g > 1 : introducing any such component will destroy 
the essential property that the evolution in f) is independent of u that 
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we require to use the normal form to project initial conditions. Thus 
we cannot use this freedom. 

• Consider (jH^ with q = 0, that is {H',T-C') = (1,0) : introducing such 
a component, multiplied by an arbitrary function of I) and e, allows 
us to widely vary the relationship between the original variable r] and 
the transformed variable i). If we maintain the exact identity between 
these variables, = f) when u = , then the model variable f) has the 
same physical meaning as the variable rj. Although not essential, this 
seems desirable as it enhances physical interpretation. We choose this 
option. 

However, the normal form transformation may instead be used to im- 
mediately put the evolution equations on the centre manifold in one of 
the well established canonical normal forms |HJ Chapt. 6], [TUl Chapt. 3] 
or [3 Chapt. 3]. This alternative also uses the freedom identified here, 
but reduces the physical connection and is not implemented here. 

• Consider with q = , that is (U',U') = (1, —1) : we cannot use this 
freedom as, by introducing terms only in f) and e into W, it destroys 
the requirement that u = is the centre manifold of the u evolution. 

• Consider (j35|) with q > 1: these homogeneous solutions allows us to 
alter the precise meaning of the transformed variable u . Note that the 
linearised dynamics about the centre manifold u = are unaffected 
by any such change as terms linear in u, corresponding to g = 1 , do 
not affect the evolution W. For example, we could use this freedom 
to replace the strict linearity condition (|28p for the evolution of u and 
instead require the normal form 

M = u + f/(t),e), and W([),0,e) = 0, (36) 

so that u measures precisely the departure of the state of the system 
from the centre manifold. For example, here the corresponding alter- 
native normal form transform is 

7] = i) + e{l + i))u + emi + i)y 
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- e^Kl + f))(12 - + 6f))u + C(e^) 



(37) 
(38) 
(39) 
(40) 



u = u-ef) + e3(3 + f^)f3 + 0(e^) , 

V = o-eu + e2f3 + C(e^), 

p = p-e^i)-e%l + [))u + 0{e^) , 



with corresponding evolution 



ut = -[1 + e\2 + [))]u - e'lil + l)^ + 0{e') , 



(41) 
(42) 
(43) 



D = p = 0. 



We adopt this ahernative when constructing the normal form of the 
dynamics of thin fluid films. 

Using such freedom we find the centre manifold and the projection of initial 
conditions is the same. For example, the direction of the projection of initial 
conditions, the tangent to the isochrons jHl E! , at the centre manifold u = 
is the same, namely 



This equivalence is reassuring because the different possible transformations 
still represent the same dynamical processes. 

3 Normal form of thin fluid fllm equations 

In this section we take the equations for a thin fluid film and transform them 
into a normal form that separates the dynamics of the viscously decaying 
modes from the large scale mode of slow evolution of the fluid film's thickness. 



"r/J re(l + [^)_e3(l + [))(2 + [)y 

Uu ^ 1 

Pu -e'(l + [)) 




(44) 
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3.1 Nondimensionalise the fluid film equations 

Consider the two-dimensional flow of a thin film of Newtonian fluid along 
a flat substrate. We adopt a nondimensionalisation based upon the char- 
acteristic thickness of the film H, and some characteristic velocity U: for 
a specific example, in a regime where surface tension drives a flow against 
viscous drag the characteristic velocity isU = a/ fi and thus the Weber num- 
ber We = a/{Ufi) = 1 ; alternatively U might be chosen to characterise the 
difference between a given initial lateral velocity and that corresponding to 
the lubrication model (^. Reverting to the general case, the reference length 
is H, the reference time H/U, and the reference pressure {lU/H. The varying 
free surface is located at y = T]{x,t), where x and y are lateral and normal 
coordinates respectively. The flow, with velocity q = [u, v) and pressure p, 
is governed by the nondimensional incompressible Navier-Stokes equations 

Re(qt + q ■ Vq) = -Vp + V^q , (45) 

where Re = UHp/fi is a Reynolds number characterising the importance 
of the inertial terms compared to viscous dissipation, supplemented by the 
continuity equation 

V-q = 0, (46) 
non-slip boundary conditions on the substrate 

q = on y = , (47) 
and tangential stress and normal stress conditions on the free surface 

(1 - T]x)iuy + v,^) = 2r]^{u^ -Vy) ony = r], (48) 

^ + (1 + ^2)3/2 = ^""y + ~ '^"^''^ + ^'"^J on y = r]. (49) 

We close the problem with the kinematic condition relating the velocity of 
the fluid on the surface to the evolution of the free surface: 

r]t = -dx f udy = v - ur]x on y = r] . (50) 
Jo 
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The fluid fllm is assumed to be so thin that the gravitational force in the 
momentum equations is neglected in this initial research project. The other 
main assumption we make is that the lateral variations are slow. 

3.2 Introduce the normal form transform 

Centre manifold techniques construct a model of the fluid flow assuming 
that lateral derivatives, dx, are small fTJ E[ e.g.]. Here we limit our aim 
to encompass flows with a lateral velocity which is also a small departure 
from the centre manifold of slow evolution. In the normal form we thus 
seek flelds u{x,y,t), \j{x,y,t), p{x,y,t) and l){x,t) which are a near identity 
transform of the physical fluid flelds. We label these "quasi-" because they 
are approximately the same as the well known physical flelds. Based upon 
the flow states we know for thin fluid fllms: 

V = i) + 0{u'/' + d^), (51) 

t; = t) + 0(u^/3 + 9^), (52) 

p = p - We f),, + O (u^/3 ^ ^ (53) 

u = u + We{[)y - |y')fi... + 0{dl) . (54) 

As in Section 1231 see equation (jSHI), we have here deflned the lateral quasi- 
velocity u to be exactly the lateral velocity fleld's departure from that of 
the centre manifold model, approximately We{l)y — \y'^)^xxx ■ Because of the 
algebraic components in the fluid fllm equations we know the normal form 
has the two trivial algebraic equations 

= and p = . (55) 

Initially we assume the normal form also has trivial algebraic boundary con- 
ditions: 

u = = Oon?/ = 0, — = p = Oony = l); (56) 
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but later we find it is necessary to modify one of tliese. Lastly, the lateral 
quasi-velocity and the fluid quasi-thickness evolve according to 

^("''^) = ^|^ + ^(^''' + ^3' ^^^^ 

nii)) = -lii)%..). + 0{d',), (58) 

where W(0, f)) = . This last property of U ensures that the trivial u = 
forms the centre manifold of the dynamics. Hence ()57|1 encapsulates that 
lateral shear flows dominantly decay by viscous diffusion, and the classic 
thin film model, f)j ~ — |(f)^f)xxx)x , is then recovered. This is our purpose 
for the normal form transformation (j51fl3^ . 

Note that the very small difference, O (u^^^ + 9^) , between the physical 
fluid thickness rj and the quasi-thickness f) identified in (j51|) means that 
evaluation on y = rj and 2/ = f) are almost everywhere interchangeable to the 
order of accuracy determined in this section. 

Also note that the form of many of the order of errors in the above are 
equivalent to having just one scaling parameter and assuming the lateral 
quasi-velocity field |u| scales with d^.'^ However, here we are primarily in- 
terested in determining the dominant nontrivial effects of a lateral velocity 
field and thus we do not seek higher order in dx except where u is involved. 
Consequently we later move to refine the approximations to a form where |u| 
is asymptotically larger and scales with — adopting a single scaling param- 
eter would hinder this development so we adopt this more flexible expression 
of asymptotic errors. 

In principle we proceed to use an iteration process to deduce the details of 
the normal form transformation. However, for the moment we only attempt 
just a little more than the first iteration. This is enough to discover highly 
nontrivial properties of the projection of initial conditions. 

I use the notation that an asymptotic error to denote the error may involve 

terms in a'"/?" for p/m + q/n > 1 . Thus, for example, here with small quantities u and dx, 
Uxx, t)^u and l)ll)xxx are all ©(u^/^ + d^). 



du 

dt 

dt 
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3.3 Fluid conservation evolves the free surface 



Seek a correction to the free surface rj 
fluid PH):^ 



i) + rj' , using conservation of 



1-^ 



di] ^ dr]' _ ^ 
dt^ dt~ " 



XI udy 





+ 0{u' + dl) 



drj' 
~dt 



Jo 



(59) 



We used tlie known evolution of the film thickness (jKH|) . that i)t ~ ~l{i)^i)xxx)x ] 
if we had not known this already, then at this step we would have discovered 
it was necessary; we can only put into rj' terms which involve u because terms 
which only involve f) are rendered ineffective in ^ by the very slow evolution 
of f). Due to the form of the right-hand side for ^ in (j59|) . now try a change 
to the fluid thickness involving a weighted integral of the lateral velocity: 



V = 

Or]' 
'dt 



'o 



{y-i))udy + 0{u' + dl) 



dt 







negligible negligible 

w^^dy\ + 0{n' + d',) 



Re dy"^ 



■^The same result is also obtained from the kinematic condition, rjt = v — urjx ony = rj , 
using the next approximation for the normal velocity field which is obtained in the next 
subsection. 
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dt " 1 Re 



du , 
w— w u 

dy 



J 



where primes on the weight function w denote y derivatives. Equate this 

dy 



w = Re{i)y - ^y^) 



to , recognising that u = on |/ = and |^ = on ?/ = f) , to deduce the 
weight function satisfies 

w" = -Re, w{0)=0, w'{i)) = 

Thus the free surface 

rj = i) + Red, [\[)y - \y^)udy + 0{v? + dl) , 
Jo 

in terms of the normal form fields f) and u. 



(60) 



(61) 



Initial conditions: Revert and (jMj) to 
rv 

f) ^ - Re 9^ / {r]y- \y'^)u dy , and u f» m - {i]y - ^y^)r]^^r, . 
Jo 

Then recall that in the normal form (fH7HHH|l . the quasi- velocity u ^ quickly 
through viscosity. Thereafter, when u = , the fluid thickness r] and the 
quasi-thickness P) are identical, by Further, r] and f) evolve according 

to the same models, (P) and (j58|) respectively, and in the normal form (|57l - 
EHI) the evolution of the quasi-thickness f) is independent of quasi-velocity. 
Hence, from the above reversion, to make the correct long term forecast with 
the lubrication model we must start with the initial fluid thickness //(x, 0) as 
claimed in (j2)). 

In the remainder of this section we demonstrate that this initial condition 
is the first nontrivial approximation in an asymptotic expansion based upon 
the normal form of thin film fluid dynamics. 



Tony Roberts, February 8, 2008 



3 Normal form of thin fluid film equations 



22 



3.4 Continuity updates the normal velocity 

Upon substituting the current approximation, u from (jK^ and v ^ d + v' 
where = from (jSSI), the continuity equation imphes 

v = t,- ru,dy-d,[ili)y'-ly^)i),,,] + 0{u' + dl). (62) 
Jo 

As necessary, for quasi-fields D = u = 0, this reduces to the well established 
normal velocity field for the centre manifold model. 



3.5 Normal momentum updates pressure 



Seek to update p ^ p -Wei)^^ + p'. Noting u^,v^ = 0{ul + 
= O(u2 + a6), consider 



and 



Re 



dv dv dv 

— — h u— h V— 

at ox oy 



+ 



dp 



+ 



.0 



dy dx'^ dy"^ 



dp' 
dy 

p' = We9.[i(f3 - + (f^ - y)u% + p" + 0{n' + d^) 



where the superscript on u^y denotes evaluation on the substrate y = 
, soon the superscript 1) will denote evaluation on the fluid surface y — 
i) , and p"{x,t) is an integration constant to be determined from the free 
surface normal stress. The normal stress condition on y = rj is equivalent to 
evaluating on ?/ = P) to as yet negligible relative error of O (u^^^ + 5^) ; thus 



p{l + r]l)+Wer],^il~W,) 



dv 
dy 



^du 



du dv 
dy dx 



+ 0{v? + dl) 



- We(f),.,. + f)f.f),.,) + / + We(l),, - if)^f) 
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= 2 [-u^. - ld^i):)X..) + - f)xu; - 0] + 0(u' + dl) 
p" = Wea,([)2{),,,.) + We |f)^(),, - 2ul - 2u;t), + 0{^x^ + ^3 , 



Hence the more refined description of tlie pressure field is 



p = p_We{t),.,.-|f)2t),, + l4[([)2 + 2f)|/-y2)[^^^^]j 
+ (f)-yK,-29.(u'^) + 0(u2 + a3. 



(63) 



3.6 Lateral momentum determines lateral velocity 

Seek to update m ~ u + m' and |j ~ + ^' ■ However, as flagged earlier, 

we are more interested in the effects of the lateral velocity field u in the 
normal form rather than higher orders in the lateral gradients dx modifying 
the leading order evolution ()58j) . Thus we change our expressions of errors to 
a form equivalent to scaling the lateral velocity u with rather than with 
as done so far. Then working to errors 0(it'^'^ + d^) enables us to resolve 
terms such as lateral diffusion u^x while omitting fifth derivatives of the free 
surface. 

Using computer algebra, see for the Reduce code, which also con- 
firms that the residual in the Navier-Stokes equations is 0(u^/^ + 9^, 
we determine 




du du duA dp d'^u d'^u 



dv! d'^v! 

=^ Re W + Re — - = Uxx - <9x[(t) - y)u°.J 
+ 2uL-2f)2.u;, + 0(u^/2 + a3, 



(64) 



At the free surface the tangential stress supplies 



f du dv\ ( dv du\ 
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^ ^= u^xdy + Au]\)^ + 0{u'/^ + dl), (65) 
Jo 

when evaluated on the free surface, y ~ f) . We could satisfy this stress 
condition by changing the velocity field in the interior through a compo- 
nent in u' . However, recall that the lateral quasi- velocity u has more direct 
meaning when the normal form is chosen so that the physical velocity u is 
independent of u except for the leading term, see (jSH). Thus, satisfy this 
inhomogeneous BC by changing the surface boundary condition (j56p for u to 

(9u f^ 

oy Jo 

Consequently a new term on the left-hand side of ()65|1 cancels the right-hand 
side terms, leaving the homogeneous BC^ = Oony = f). Now consider 
the lateral momentum update equation (|64p : all terms in the right-hand 
side involve u and so they are all placed into the evolution W . Thus the 
lateral velocity is still given by (jSH), but we improve the description of the 
evolution to 

Re ^ = ^yy + - ^^[(f) - vXy] + 2uL - 2f)^u;^ + 0{xV'l^ + dl) . (67) 

As required, the normal form equation ()67|1 with boundary condition ()66|1 
has u = as an attractive invariant (centre) manifold. Then with lateral 
quasi-velocity u = the lateral velocity field reduces to the conventional 
thin fluid film approximation u ~ We(f)?/ — \y'^)^xxx ■ 



3.7 Fluid conservation refines the evolution 

We seek a further refinement to the description of the fluid thickness 1)6111 
using conservation of fluid (jSUJ. The aim is to discover more effects of the 
lateral velocity and so we work to errors 0{y?+d^ . There is considerable de- 
tail in determining the new terms, relegated to the computer algebra in [TB] , 
but the basic technique follows that used in Section KlHI To the required 
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order of accuracy we know the quasi-thickness evolution ()58|1 : thus the only 
freedom available is to update the fluid thickness from by some small 
change rj". Substitute the lateral velocity field into to obtain (using 
superscripts and f) to denote evaluation on the bed and quasi-surface, but 
now also using superscript y to denote ■ dy) 

^S' ^xxxy ' qV vx^xxy ' 3 V ')xx^xy ' 'J Vx^xy 3 'J ^xxx 

+ 4l)'i)lnly - 2f}().(uL)^ - ml4 + «fj^ + [).(ufj^ 

+ 0{u^ + d^,). (68) 

All these residual terms are (9(m(9^), there is no component of 0(u(9^). We 
summarise the details of the derivation in the next paragraph, but a solution 
for the above refinement of the fluid thickness rj is 

V'(f3,u) = Re |/J(^y^-l[)y3_ 1(32^2^ 1^3^ _|[^4)^_^^ 

+ i)4li-ly' - W + Wv - l^^xxdy 

+ WJli^y - f))u.. dy + j;(2y - if))u,. dy 

+ m^\)iui + f fi^^f^.fi^.u'' + i\)^\)iv^ + ?,\)%Ax 

+ 2f)^[)..u^ + [)^()...u^ - 2\)%luly] . (69) 

Thus for this rj" the normal form for the fluid thickness is now 

r/ = () + Re9j;(t)y - \y^)udy + V'(f), u) + ^(u^ + dl) . (70) 

In the next section we use this transformation to determine more details 
about the initial conditions for the model of thin fluid film (P). 

But before proceeding we overview the machinations needed to derive 
the refinement (j69p . Recall that the left-hand side ^ is approximated by a 
formal expression: 

drj' 1 drj' d'^u 
dt Re du dy"^ ' 
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Thus the components in the residual on the right-hand side of fl68|l which 
involve two or more derivatives of y can be matched by a component in the 
refinement rj" with two less derivatives in y. The other components, involving 
the first ^/-derivative or integrals of the quasi- velocity u, must be matched by 
an integral component in the refinement f]" as we did in Section I achieve 
such matching through considerable trial and error, and in three stages. 

• Starting with the term with the most y-integrals, namely (u^^^.)'', the 
three ^/-derivatives in this boundary contribution is matched through 
multiple integration by parts of dfUxx oc u^xyy with a weight function 
which is quartic in y. This and the bed component in f)^u°^^y gives the 
first integral in (j69p . This first term in (j69|) also changes some other 
terms in the residual (|68|) . Then we progress to the other terms with 
less ?/- integrals and match them with simpler integrals as seen in ()69j) . 

• When all components as yet unmatched have two or more y-derivatives 
then we simply directly match them with the boundary contributions 
seen in fjHTH) . 

• However, there is a difficulty: terms with odd y derivatives and eval- 
uated at the quasi-surface cannot be matched in the above procedure. 
Instead we use the surface boundary condition (jHS)), either directly or its 
derivatives, to eliminate such terms. For example, differentiating (j66p 
with respect to x gives 

= Jluxx. dy + 5i)xul + 4i)xA + 4f)'u^, + O [n'/' + d^) 
= Jluxx. dy + 5i)xul + Ai)xxul - + O (u^/^ + ^ 

where the last term u^y has a single y-derivative and so is replaced by 
its leading order approximation, namely —i)xU^y. The above and higher 
order derivatives are implemented in computer algebra |13^. Similarly, 
time derivatives of (jUnj) . recalling Reuj = Uyy , provide replacements for 
higher order odd y-derivatives of u. 
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With these methods we find that is the solution for the update 77". 

4 Revert the normal form to give initial con- 
ditions 

Given an initial state of the fiuid film we determine the appropriate initial 
condition for the fiuid thickness in the lubrication model (^. First translate 
the initial condition to the normal form variables f) and u, then since the 
evolution of f) is independent of quasi- velocity u, the correct initial condition 
for the lubrication model ((H) is ri{x, 0) = f)o • 

Suppose the fiuid has initial thickness rjo^x), lateral velocity field uo{x, y), 
normal velocity field vo{x, y) and pressure field Pq{x, y). The fields fo and po 
must be consistent with the fiuid equations ()45flKn|l . but then play no role in 
the following. We revert (j54j) to obtain the initial quasi-lateral velocity 

uo = uo - We(r7oy - \y^)mxxx + 0{v? + dl) , (71) 

as the quasi-thickness f) = 77 + Oiy?^"^ + 9^), by (fTOI) for example, and so 
f) is replaced by the initial fiuid thickness rjQ to our order of accuracy. This 
quasi-lateral velocity Uq characterises the distance from the initial fiuid state 
to a state of the lubrication model. 

Now find how this distance affects the initial fluid thickness. Revert ()7()|l 
and evaluate at the initial fluid state, using (|7T|l and f) = + 0{y?/'^ + (9^) , 
to give the initial quasi-thickness 

f)o = r/o - Re9, / (r/oy - \y^)uody - v"{Vo,no) + 0{u' + 9^) , (72) 
Jo 

where 1]" is specifled by Recall the purpose of the normal form is to 

ensure the evolution of quasi-thickness f) is independent of the quasi- velocity 
and hence the transient viscous decay of lateral velocity will bring the fluid 
to the solution of the lubrication model (P) which started from the initial 
condition ri{x, 0) = P)o as specifled in (f72|) . This initial condition permits us 
to make accurate long term forecasts with the model. 
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5 Conclusion 

We use a normal form transformation to illuminate the dynamics of a thin 
layer of fluid. This is achieved by decoupling the slow long-lasting lubri- 
cation mode from the viscously decaying lateral shear modes. A simple 
example shows that the differential algebraic nature of the fluid equations 
is handled by straightforward modifications of the usual procedure to con- 
struct a normal form. Further, invoking the slowly-varying assumption that 
lateral derivatives are small enables us to deduce a normal form for the fluid 
dynamics albeit with significant technical detail requiring computer algebra 
to check. This normal form then provides us with the rationale to choose the 
initial condition Q to make forecasts with the lubrication model(P]). The 
next challenge is to connect this normal form analysis to the projection of 
initial conditions, started in jTHj, which is based upon the adjoint near a 
centre manifold model. 

Acknowledgement: I thank Sergey Suslov for his valuable input into all 
stages of this work. 
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